Cell Reports Methods

Fluorescent in situ hybridization (FISH) enables highly sensitive, high-resolution detection of gene transcripts. Moreover, by employing multiple probes, this technique allows for multiplexed, simultaneous detection of distinct gene expression patterns spatiotemporally, making it a valuable spatial transcriptomics approach. Owing to these advantages, FISH techniques are rapidly being adopted across diverse areas of basic biology. However, conventional protocols often rely on volatile, toxic reagents such as formalin or methanol, posing potential health risks to researchers. Here, we present a safer protocol that replaces these chemicals with low-toxicity alternatives, without compromising the high detection sensitivity of FISH. We validated this protocol using both in situ hybridization chain reaction (HCR) and signal amplification by exchange reaction (SABER)-FISH in frozen sections of various model organisms, including mouse (Mus musculus), amphibians (Xenopus laevis and Pleurodeles waltl), and medaka (Oryzias latipes). Our results demonstrate successful multiplexed detection of morphogenetic and cell-type marker genes in these model animals using this safer protocol. The protocol has the additional advantage of requiring no proteolytic enzyme treatment, thus preserving tissue integrity. Furthermore, we show that this protocol is fully compatible with EGFP immunostaining, allowing for the simultaneous detection of mRNAs and reporter proteins in transgenic animals. This protocol retains the benefits of highly sensitive, multiplexed, and multimodal detection afforded by integrating in situ HCR and SABER-FISH with immunohistochemistry, while providing a safer option for researchers, thereby offering a valuable tool for basic biology.

Circulating tumor cells (CTCs), and especially CTC-clusters, are linked to poor prognosis and may reveal mechanisms of metastasis and treatment resistance. Therefore, developing unbiased methods for the functional characterization of CTCs in liquid biopsies is an urgent need. Here, we present an evaluation of multiplex imaging mass cytometry (IMC) to analyze CTCs in mice with human xenograft tumors. In a single-step process, IMC uses metal-labeled antibodies to simultaneously detect a large number of proteins/modifications within minimally manipulated small volumes of blood from the tail vein or heart. We used breast cancer cell lines and a patient-derived xenograft (PDX) to assess antibodies for cross-species interpretation. Along with manual verification, HALO-AI-based cell segmentation was used to identify CTCs and quantify markers. Despite some limitations regarding human-specificity, this technology can be used to investigate the effect of genetic and pharmacological interventions on the properties of single and cluster CTCs in tumor-bearing mice.

Background: The enterotype concept proposed that gut microbiomes cluster into discrete types, but subsequent critiques demonstrated that such clustering depends on methodological choices, that the number of clusters is not fixed, and that faecal samples cannot capture spatial heterogeneity along the gastrointestinal tract. The stomach remains particularly understudied, and no systematic classification exists for gastric microbial community types. Methods: We assembled a multi-cohort dataset of 566 gastric mucosal samples spanning healthy controls to gastric cancer, with both Helicobacter pylori (HP)-negative and HP-positive individuals. Critically, we applied the key methodological lessons of the enterotype debate: we used a variational autoencoder (VAE) for dimensionality reduction to learn a continuous latent representation without forcing discrete structure, determined the optimal number of clusters using the Silhouette index (an absolute validation measure) across K=2 to K=10 rather than arbitrarily selecting a cluster number, and performed transparent evaluation of multiple clustering solutions. This VAE-plus-silhouette workflow directly addresses the critiques leveled against the original enterotype analysis. Results: Four gastotypes were identified, with K=4 achieving the highest mean silhouette score, indicating good cluster cohesion and separation. Two gastotypes (Variovorax-type and Trabulsiella-type) were significantly enriched in HP-positive samples, while two gastotypes (Bacteroides-type and Streptococcus-type) were significantly enriched in HP-negative samples. Random Forest and Gradient Boosting achieved excellent baseline performance for predicting HP infection (AUC = 0.990 and 0.993). Conclusions: The VAE-plus-silhouette workflow provides a robust, data-driven approach for identifying gastotypes without forcing discrete structure or arbitrarily fixing cluster numbers. Using this framework, we identified four gastotypes with significantly different HP infection rates. Variovorax-type and Trabulsiella-type showed strong HP-positive enrichment, while Bacteroides-type and Streptococcus-type showed strong HP-negative enrichment. These findings demonstrate that methodological advances from the enterotype controversy can be successfully transferred to the stomach, offering a reproducible taxonomy for stratifying HP infection status with potential clinical utility.

Background. Pancreatic ductal adenocarcinoma (PDAC) has a five-year survival rate of approximately 12%, largely because it is typically diagnosed at an advanced stage. CT-based computational methods for early detection exist but rely on black-box deep learning or large texture feature sets without tissue-specific interpretability. Methods. We developed Virtual Spectral Decomposition (VSD), which applies six parameterized sigmoid functions S(HU) = 1/(1+exp(-alpha x (HU - mu))) to standard portal-venous CT, decomposing each pixel into tissue-specific response channels for fat (mu=-60), fluid (mu=10), parenchyma (mu=45), stroma (mu=75), vascular (mu=130), and calcification (mu=250). Dendritic Binary Gating identifies structural content per channel using morphological filtering, enabling co-firing analysis and lone firer identification. A 25-feature signature was extracted per patient. Three independent datasets were analyzed: NIH Pancreas-CT (n=78 healthy), Medical Segmentation Decathlon Task07 (n=281 PDAC, paired tumor/adjacent tissue), and CPTAC-PDA from The Cancer Imaging Archive (n=82, multi-institutional, with DICOM time point tags). The same six sigmoid parameters were used across all datasets without retraining. Results. VSD achieved AUC 0.943 for field effect detection (healthy vs cancer-adjacent parenchyma) and AUC 0.931 for patient-stratified tumor specification on MSD. On CPTAC-PDA, VSD achieved AUC 0.961 (6 features) and 0.979 (25 features) for distinguishing healthy from cancer-bearing pancreas on scans obtained prior to pathological diagnosis. All significant features replicated across datasets in the same direction: z_fat (d=-2.10, p=3.5e-27), z_fluid (d=-2.76, p=2.4e-38), fire_fat (d=+2.18, p=1.2e-28). Critically, VSD severity did not correlate with days-from-diagnosis (r=-0.008, p=0.944) across a range of day -1394 to day +249. Patient C3N-01375, scanned 3.8 years before pathological diagnosis, had VSD severity 1.87, well above the healthy mean of 0.94 +/- 0.33. The tissue transformation signature was temporally stable, indicating an early, persistent tissue state rather than a progressively worsening process. Conclusions. VSD with Dendritic Binary Gating detects a stable pancreatic tissue composition signature on standard CT that is present years before clinical diagnosis, validated across three independent datasets without parameter adjustment. The six sigmoid channels map to biologically meaningful tissue components through a fully transparent interpretability chain. The temporal stability of the signal implies a detection window of 3-7 years, consistent with known PanIN-3 microenvironment transformation timelines. VSD functions as a single-scan screening tool applicable to any abdominal CT performed during the pre-clinical window.

Lightweight epidemic calculators are widely used for teaching and rapid scenario exploration, yet many omit the methodological detail needed for scientific reuse. We present a browser-native SIR calculator that exposes forward Euler and classical fourth-order Runge--Kutta (RK4) integration alongside epidemiologically interpretable outputs and a population-conservation diagnostic. The implementation is anchored to analytical properties of the deterministic SIR system, including the epidemic threshold, the peak condition, and the final-size relation. Benchmark experiments show that RK4 is essentially step-size invariant over practical discretizations, whereas Euler at a coarse one-day step overestimates peak prevalence by 3.97% and final size by 0.66% relative to a fine-step RK4 reference. These results demonstrate that browser-based tools can support publication-quality computational narratives when solver choice, diagnostics, and assumptions are treated as first-class outputs.

PurposeAntibodies to Epstein-Barr virus (EBV) proteins can predict nasopharyngeal carcinoma (NPC) risk. We previously defined a prototype EBNA1 protein panel and multiplex immunoblot assay that distinguishes NPC risk several years pre-diagnosis. Assay throughput and specificity are critical to effectively implement a population-level screening program. Here, we developed a strip test assay - EBNA1 SeroStrip-HT - with an objective to increase throughput and maximize specificity. Experimental DesignEBNA1 full-length (FL) and glycine-alanine repeat deletion mutants (dGAr) were purified from insect and mammalian cells to screen serum IgA/IgG from prospective cohorts in Singapore and Shanghai, China, with known time intervals to NPC diagnosis. Twenty pre-diagnostic sera within 4 years to diagnosis were compared to 96 healthy controls using a nested case-control study design. ResultsIgA to mammalian-derived EBNA1 dGAr achieved 85.0% sensitivity and 94.8% specificity (AUC, 0.939) for NPC status. IgA to insect-derived EBNA1 dGAr showed the same sensitivity (85.0%) and similar specificity (93.8%) (AUC, 0.941). IgA to insect-derived EBNA1 FL had a higher 90% sensitivity, but lower 91.7% specificity (AUC, 0.940). Combining EBNA1 FL and dGAr results showed that subjects positive for both proteins had a 243.67 odds ratio for NPC incidence compared to double-negative scores. ConclusionThis study demonstrated the efficacy of EBNA1 SeroStrip-HT for NPC risk assessment and stratification in high- and intermediate-risk populations, yielding high accuracy and a 12-fold increased throughput over the prototype. The insect system was appropriate for large-scale production of purified EBNA1. Larger, geographically diverse cohorts are warranted to confirm these results, especially in low-incidence populations.

End-stage liver disease (ESLD) is one of the leading causes of death worldwide. Currently, the only curative option for patients with ESLD is liver transplantation. However, the demand for donor livers far exceeds the available supply, partly because many potentially viable livers are discarded following biopsy evaluation. While biopsy is the gold standard for assessing liver histological features related to graft quality and transplant suitability, it often leads to high discard rates due to its susceptibility to sampling errors and limited spatial coverage. Besides, biopsy is invasive, time-consuming, and unavailable in clinical facilities with limited resources. Here, we present an AI-assisted photoacoustic/ultrasound (PA/US) imaging framework for quantitative assessment of human donor liver graft quality and transplant suitablity at the whole-organ scale. With multimodal volumetric PA/US images as the input, our deep-learning (DL) model accurately predicted the risk level of fibrosis and steatosis, which indicate the graft quality and transplant suitability, when comparing with true pathological scores. DL also identified the imaging modes (PAI wavelength and B-mode USI) that correlated the most with prediction accuracy, without relying on ill-posed spectral unmixing. Our method was evaluated in six discarded human donor livers comprising sixty spatially matched regions of interest. Our study will pave the way for a new standard of care in organ graft quality and transplant suitability that is fast, noninvasive, and spatially thorough to prevent unnecessary organ discards in liver transplantation.

ELF4 is an ETS family transcription factor involved in immune regulation, and germline loss-of-function mutations in ELF4 have been known as deficiency in ELF4, X-linked (DEX). To date, ELF4-related disease has been exclusively associated with germline mutations. Here, we report a pediatric patient with recurrent mucocutaneous inflammation and periodic fever caused by a somatic truncating mutation in ELF4. By directly comparing ELF4-mutant and wild-type immune cells within the same individual using full-length single-cell RNA sequencing, we identified mutation-associated transcriptional alterations across multiple immune cell types. Pathway analyses revealed cell type-specific immune alterations, characterized by reduced antiviral and interferon-related signaling in NK cells and enhanced inflammatory pathways related to Th17 differentiation and inflammatory bowel disease in CD16 monocytes. This study expands the disease spectrum of ELF4 deficiency by identifying somatic truncation of ELF4 as a genetic mechanism underlying autoinflammatory diseases and biased immune programs.

Ewing sarcoma (EwS) is an aggressive, human-exclusive tumor typically driven by the EWS::FLI1 fusion protein. To assess whether the neomorphic functions of EWS::FLI1 are fundamentally dependent on evolutionarily recent cofactors such as ETS transcription factors (ETS-TFs), Plycomb group (PcG) proteins, CBP/p300, or specific subunits of the BAF complex, we expressed EWS::FLI1 in the model organism Saccharomyces cerevisiae. This minimal system was chosen because several key EWS::FLI 's cofactors possess greatly reduced sequence homology (e.g., BAF) or are lacking altogether (e.g., ETS-TFs, PcG, or CBP/p300). We used co-IP/MS to map the yeast interactome, Chip-Seq to identify gDNA binding sequences, RNA-Seq for global gene expression, and engineered reporters to test conversion of (GGAA) tandem repeats (GGAASat) into neoenhancers. We found that the yeast EWS::FLI1 interactome was more limited and qualitatively distinct from its human counterpart, sharing core machinery (e.g. RNA Polymerase II, FACT) but lacking the BAF/SWI-SNF and spliceosome complexes, and showing strong enrichment for the SAGA chromatin remodeling complex. We also found that EWS::FLI1 binds to hundreds of sites in the yeast genome with a clear preference for putative ETS-TF consensus sequences and (CA) dinucleotide repeats. Yet, EWS::FLI1 expressing cells presented only minimal transcriptional dysregulation, a stark contrast to the extensive changes observed in humans and Drosophila cells. Finally, we found that EWS::FLI1 successfully converted silent GGAASat sequences into active enhancers in yeast. This remarkable result occurs despite the absence of homologs for key human activators, such as CBP/p300, strongly suggesting that EWS::FLI1 can mobilize functionally related, non-homologous pathways to establish neoenhancers at GGAASat sites. Altogether, our results indicate that EWS::FLI1's core ability to drive GGAASat-dependent gene expression is a conserved, ancient property, while GGAASat-independent extensive transcriptome reprogramming is dependent on co-factors and pathways specific to animal cells.

Streptococcus pneumoniae is the leading cause of empyema and pneumonia in children, and monitoring of effectiveness of polyvalent pneumococcal vaccines has been essential for controlling invasive pneumococcal disease (IPD) in children and elderly adults. Conventional serotyping of pneumococci has relied on Quellung reaction following laboratory culture, however more recently whole genome sequencing (WGS) has been implemented in many reference laboratories to enhance traditional typing. Pleural fluid samples from cases with empyema are often culture negative, limiting the utility of WGS and requiring polymerase chain reaction (PCR) or 16S rRNA sequencing to detect S. pneumoniae. These molecular methods have limited sensitivity and capacity to characterise pneumococcus in clinical samples, especially in specimens with a low pathogen abundance. This study applied capture-based enrichment (tNGS) to identify and characterise S. pneumoniae directly from pleural fluid samples. A total of 51 pleural fluid samples were subjected to tNGS with a custom probe panel, for 39 known positive fluids collected from IPD cases between 2018-2025 in New South Wales, Australia. tNGS results were benchmarked against molecular-based serotyping. Our tNGS achieved 100% sensitivity and specificity in detecting S. pneumoniae. Serotyping results were concordant with PCR and 95% (37/39) of S. pneumoniae PCR positive pleural fluid cases could be serotyped using tNGS. Standard molecular methods however could only determine serotype in 56% (22/39) of samples. This tNGS enabled 39% improvement in ability to directly identify and serotype IPD-associated serotypes of S. pneumoniae in difficult-to-culture pleural fluids can significantly enhance laboratory surveillance of IPD as well as our understanding of vaccine effectiveness.

Obtaining novel antibodies against specific protein targets is a widely important yet experimentally laborious process. Meanwhile, computational methods for antibody design have been limited by low success rates that currently require resource-intensive screening. Here, we introduce Germinal, a broadly enabling generative pipeline that designs antibodies against specific epitopes with nanomolar binding affinities while requiring only low-n experimental testing. Our method co-optimizes antibody structure and sequence by integrating a structure predictor with an antibody-specific protein language model to perform de novo design of functional complementarity-determining regions (CDRs) onto a user-specified structural framework. When tested against four diverse protein targets, Germinal successfully designed functional antibodies across all targets and binder formats, testing only 43-101 designs for each antigen. Validated designs also exhibited robust expression in mammalian cells and high sequence and structural novelty. We provide open-source code and full computational and experimental protocols to facilitate wide adoption. Germinal represents a milestone in efficient, epitope-targeted de novo antibody design, with notable implications for the development of molecular tools and therapeutics.

Effective predictive biomarkers for immune checkpoint inhibitor (ICI) therapy remain an unmet need across solid tumors. Here, we present an integrated spatial proteomics workflow that combines in situ proximity ligation assay with multiplexed immunofluorescence to directly resolve PD1/PDL1 signaling events at the level of defined cellular phenotypes and their spatial organization within intact tumor tissue. Applied as a proof of concept to tumor samples from patients with metastatic urothelial carcinoma treated with pembrolizumab, this approach reveals that PD1/PDL1 interactions specifically involving cytotoxic CD8CD3 T cells are significantly enriched in complete responders, while such interactions are rare in patients with progressive disease. This interaction defined T cell subset achieves superior discrimination of clinical response compared to single marker PDL1 expression or immune cell abundance alone. By integrating direct detection of protein protein interactions with high dimensional single cell phenotyping, our workflow provides a mechanistically informed, spatially resolved biomarker of functional immune engagement. Beyond urothelial carcinoma, this platform establishes a generalizable framework for translating spatial signaling biology into predictive tools for immunotherapy response across tumor types.

Mycetoma is a neglected tropical disease caused by various bacterial and fungal pathogens that has a significant health impact across a broad geographically defined "mycetoma belt" spanning South America, Africa and Asia. Histologically, mycetoma is characterised by invasive and destructive granuloma development in the skin, deep tissues and bone, leading to tissue destruction, deformities and high morbidity. The presence of macroscopic, highly compacted pathogen microcolonies, or "grains," is a key diagnostic feature, and the formation of grains supports pathogen persistence and disease chronicity. However, there is a paucity of information on immune responses in mycetoma patients and on the relative importance of phylogeny and/or grains in establishing the local immune landscape. Here, we used spatial proteomics to examine the distribution of 43 immune-related proteins in surgical biopsies from 11 patients with mycetoma of bacterial (Actinomycetoma; Actinomadura pelletierii and Streptomyces somaliensis; n=6) and fungal (Eumycetoma; Madurella mycetomatis; n=5) origin. Using mixed-effects modelling, an exploratory analysis across species and pathogen classes revealed few significant differences in immune marker expression. In contrast, and independently of pathogen class, the cellular infiltrate closest to grain boundaries had higher per-cell expression of CD66b+, ARG1, and VISTA. The preferential accumulation of CD66b+ARG1+VISTA+ cells at grain boundaries was confirmed by quantitative immunofluorescence analysis. Hence, the local tissue microenvironment surrounding the mycetoma grain represents a specialised immunosuppressive niche, with parallels to the tumour microenvironment.

Neural excitability fluctuates with sensory events, creating windows of opportunity to enhance brain stimulation. Repetitive transcranial magnetic stimulation (TMS), including intermittent theta burst stimulation (iTBS), is a promising treatment for neurological and psychiatric disorders, but does not account for fluctuations in neural excitability, likely contributing to variable outcomes. Sensory Entrained TMS (seTMS) leverages sensorimotor oscillations to enhance corticospinal responses, but the sustained effects as a repetitive protocol are unknown. We extend seTMS to iTBS, measuring motor-evoked potentials (MEPs) as a physiological readout. In a randomized crossover study comparing standard iTBS with sensory entrained iTBS (se-iTBS; n=20), we found that se-iTBS more than doubled the MEP effect (55% vs 26% MEP enhancement) and persisted for at least 30 minutes. Notably, at least 80% of participants showed larger responses with se-iTBS at all time points. se-iTBS may provide a robust and practical framework for optimizing TMS that bridges electrophysiological mechanisms and clinical applications.

Patients who participate in intracranial neuroscience research make invaluable contributions to our understanding of the brain, accelerating the development of neurotechnological interventions. Engagement of patients as part of this research presents unique challenges, where study goals can be distant from immediate clinical applications and require specialized domain knowledge. Yet methods for meaningfully integrating patient communities as part of these research efforts is essential, as intracranial neuroscience guides the application of artificial intelligence for understanding and enhancing human cognition. In order to identify what patients consider meaningful research engagement we interviewed individuals who participated in a study during their Deep Brain Stimulation (DBS) surgery and attended a group event where they interacted with our research team. Analysis of semi-structured interviews identified four main themes: interest in science and the future of clinical care, contributing to science to improve lives, connecting with others, and accessibility considerations. Based on these insights, we propose strategies for transformational participation of patient communities in intracranial neuroscience research with respect to engagement objectives, communication and scope. This approach offers a foundation for sustaining relationships between scientists and communities rooted in trust and transparency, to ensure that impacts of neurotechnology on human health and cognition are aligned with patient needs as well as desired public values.

Gene networks (GNs) encode diverse molecular relationships and are central to interpreting cellular function and disease. The heterogeneity of interaction types has led to computational methods specialized for particular network contexts. Large language models (LLMs) offer a unified, language-based formulation of GN inference by leveraging biological knowledge from large-scale text corpora, yet their effectiveness remains sensitive to prompt design. Here, we introduce Gene-Relation Adaptive Soft Prompt (GRASP), a parameter-efficient and trainable framework that conditions inference on each gene pair through only three virtual tokens. Using factorized gene-specific and relation-aware components, GRASP learns to map each pair's biological context into compact soft prompts that combine pair-specific signals with shared interaction patterns. Across diverse GN inference tasks, GRASP consistently outperforms alternative prompting strategies. It also shows a stronger ability to recover unannotated interactions from synthetic negative sets, suggesting its capacity to identify biologically meaningful relationships beyond existing databases. Together, these results establish GRASP as a scalable and generalizable prompting framework for LLM-based GN inference.

An efficacious blood-stage malaria vaccine would serve as a highly useful public health tool alongside licensed vaccines targeting the pre-erythrocytic life cycle stage of the Plasmodium falciparum parasite. RH5 is the leading blood-stage malaria vaccine candidate antigen due to its highly-conserved sequence and non-redundant role in merozoite invasion of red blood cells. Following encouraging immunogenicity data in UK and Tanzanian Phase Ia/b vaccine trials, RH5-based vaccines have progressed to Phase IIb evaluation in Burkina Faso in recent years. Here, we report a Phase Ia clinical trial in malaria-naive UK adults to assess the safety and immunogenicity of the malaria vaccine candidate RH5.1 soluble protein with Matrix-M adjuvant using two different booster dosing regimens: 10-10-10 micrograms versus 50-50-10 micrograms RH5.1, both delivered in a 0-1-6-month schedule with 50 micrograms Matrix-M adjuvant per dose (ClinicalTrials.gov NCT06141057). A total of n=24 participants were recruited to this study, with n=23 completing all follow-up visits through to 1 year following final vaccination. The RH5.1/Matrix-M formulation was well-tolerated in this population, with injection site pain, myalgia and fatigue being the most commonly reported symptoms up to 7 days post-vaccination. There were no serious adverse events, adverse events of special interest, or suspected unexpected serious adverse reactions reported over the course of the trial. Both vaccination regimens were similarly immunogenic; no differences were observed in peak anti-RH5.1 serum IgG concentrations, in vitro functional anti-parasitic activity, avidity, or durability. Our findings build on other observations from clinical trials of adjuvanted RH5.1 indicating that humoral immunogenicity can be enhanced by delaying the final booster vaccination, but that there is limited impact of fractionation of the final dose. These insights can help to guide the next steps of multi-antigen, multi-stage malaria vaccine development in malaria-endemic settings.

How embryonic cells generate large clones of cells in the adult represents a fundamental question in biology. Here, using melanocyte stem cells (McSCs) in the zebrafish as a model, we explore the function of the master melanocyte transcription factor (MITF) in safeguarding McSCs in embryonic development and their potential to pigment large clones in the adult. MITF is well known is for its role in the specification of melanoblasts from the neural crest (NC) and their differentiation into melanocytes, yet little is known about how this activity shapes the stem cell lineages. Here, we use live imaging coupled with single-cell transcriptomics and lineage tracing to show that MITF (mitfa in zebrafish) protects the melanocyte stem cell (McSC) fate in zebrafish. Utilizing a temperature sensitive mitfavc7 mutant, we show loss of Mitfa leads to a surprising premature and aberrant expansion of McSC progeny at the niche during embryogenesis, coupled with novel emergent transcriptional cell states. Linage tracing of McSCs from the embryonic to juvenile stages reveals Mitfa activity is subsequently required in regeneration by Schwann cell-like and melanocyte stem cell progenitors that serve as a reservoir for fast-responding pigment progenitors. Thus, the impact of Mitfa loss on the melanocyte lineage is cell-state and stage-specific. The emergent cell states upon mitfa loss may have important implications for our understanding the loss of MITF activity in human genetic disease and melanoma.

Detailed agent-based simulations are increasingly used to support policy decisions, but their computational cost and complex uncertainty structure make systematic scenario analysis challenging. We present a data-driven, uncertainty-aware decision support (DDUADS) workflow for using stochastic simulation models as decision-support tools under limited computational budgets. The approach combines several established techniques-sensitivity screening, Bayesian calibration using simulation-based inference, and multi-surrogate model integration for translational efficiency-into a coherent pipeline that enables uncertainty-aware policy analysis. Rather than producing a single baseline, the calibration stage yields a posterior distribution over plausible model parameterizations, allowing flexible, uncertainty-aware forward projections. We demonstrate the DDUADS workflow on the INFORM-HIV agent-based model of HIV transmission in Chicago to evaluate potential disruptions in antiretroviral therapy (ART) and pre-exposure prophylaxis (PrEP) use. While the specific application is HIV modeling, the challenges and techniques described here arise in other simulation studies and can be applied to decision support in other domains.

Cervical spinal cord injury (SCI) causes profound and persistent loss of hand function, and effective neuromodulation strategies remain limited. We report the first-in-human implantation of a 32-contact cervical epidural paddle array in two individuals with severe chronic SCI. Individualized motor pool recruitment maps, derived from systematic bipolar and multipolar configurations, enabled person-specific stimulation parameters. Optimized stimulation restored volitional hand opening, closing and coordinated upper-limb movements that were previously unattainable. This approach achieved a >91% success rate in complex reach-grasp-lift-release sequences, supported by substantial gains in range of motion, grip, and pinch strength. Electrophysiological and kinematic analyses demonstrated parameter-dependent, selective recruitment of flexor and extensor motor pools. Personalized stimulation programs integrated with goal-directed activities enabled functional hand use in home and community settings, sustained over several months of continued autonomous use. These findings establish a mechanistically grounded and translational framework for restoring upper-limb function after chronic severe SCI.